Unconventional superconducting pairing by conventional phonons 
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The common wisdom that the phonon mechanism of electron pairing in the weak-coupling 
Bardeen-Cooper-Schrieffer (BCS) superconductors leads to conventional s-wave Cooper pairs is re- 
vised. An inevitable anisotropy of sound velocity in crystals makes the phonon-mediated attraction 
of electrons non-local in space providing unconventional Cooper pairs with a nonzero orbital momen- 
tum in a wide range of electron densities. As a result of this anisotropy quasi-two dimensional charge 
carriers undergo a quantum phase transition from an unconventional d-wave superconducting state 
to a conventional s-wave superconductor with more carriers per unit cell. In the opposite strong- 
coupling regime rotational symmetry breaking appears as a result of a reduced Coulomb repulsion 
between unconventional bipolarons dismissing thereby some constraints on unconventional pairing 
in the Bose-Einstein condensation (BEC) limit. The conventional phonons, and not superexchange, 
are shown to be responsible for the d-wave symmetry of cuprate superconductors, where the on-site 
Coulomb repulsion is large. 
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A well-known quantum mechanics theorem |l| states 
that the coordinate wave function of two particles does 
not become zero (or has no nodes) in the ground state. If 
it has no nodes the wave function must certainly be sym- 
metrical with respect to the interchange of the particle 
coordinates, i"i <=^ V2] for, if it were to be antisymmetri- 
cal, it would vanish at ri = T2- Thus the system of two 
identical fcrmions must have a spin-singlet ground state, 
as the complete wave function (which includes spin coor- 
dinates) must be antisymmetrical. 

It has been thought for a long while that Cooper 
pairs in the BCS theory with the conventional electron- 
acoustic) phonon interaction (EPI) are singlets and their 
wave function is isotropic (s-wave), in apparent agree- 
ment with the theorem (2j. The interaction has been 
thought to be local in space, so it could not lead to a 
higher angular momentum pairing. If the pairing interac- 
tion is non-local, the BCS theory admits unconventional 
Cooper pairs with nonzero orbital momentum [3[ , as the 
theorem is not applied to more than two particles. Thus 
it has gone unquestioned that the unconventional pair- 
ing requires unconventional electron-phonon interactions 
with specific phonons and poor screening 0, H, H, 0> H] , 
sometimes combined with anti-ferromagnetic fluctua- 
tions 9] and vertex corrections [Toj . or non-phononic 
mechanisms of pairing (e.g. superexchange [ll|), and a 
specific shape of the Fermi surface. Since the rotational 
symmetry breaking is a many-body effect, any supercon- 
ductor should seem to be s-wave in the strong-coupling 
limit (T3 |. where pairs are individual (e.g. bipolarons 
[HI]) rather than overlapping Cooper pairs. 

Here the symmetry of the superconducting state me- 
diated by conventional acoustic phonons is reconsidered. 
The sound speed anisotropy leads to a double surprise: 
a) the BCS state of layered crystals is d-wave in a wide 
range of carrier densities, b) the strong-coupling BEC 
state can break the rotational symmetry as well as the 



weak-coupling BCS state. 

In the framework of the BCS theory the symmetry of 
the order parameter A(k) and the critical temperature, 
Tr, are found by solving the linearised "master" equation 

a. 
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The interaction V(k, k') comprises the attraction, 
—V p h(<l), mediated by acoustic phonons, and the 
Coulomb repulsion, V c (q), as 

F(k,k') = -V ph (q)Q(cJD - |&|)e(w D - l&'l) 

+ v c (q)eH-|6<|)e( Wp -|6c|), (2) 

where V^j,(q) = C 2 /NMcf is the square of the matrix el- 
ement of the electron-phonon interaction divided by 
the square of the acoustic phonon frequency, o-> q = ciq, 
ci is sound velocity, M is the ion mass, N is the num- 
ber of unit cells in the crystal, and £k is the electron 
energy relative to the Fermi energy. The deformation 
potential matrix element C is roughly the electron band- 
width in rigid metallic or semiconducting [l5| lat- 
tices, q — k — k' is the electron momentum transfer with 
the magnitude q = 2 1 / 2 fc,p[l — cos^] 1 / 2 , where $ is the 
angle between k and k', and hkp is the Fermi momen- 
tum. Thcta functions in Eq.([T]) (Q(x) — 1 for positive x 
and zero otherwise) account for a difference in frequency 
scales of the electron-phonon interaction, wrj, and the 
Coulomb repulsion, uj p 3> lod, where cud and uj p are the 
maximum phonon and plasmon energies, respectively. 

If one neglects anisotropic effects, replacing V p h{(\} and 
V c (q) by their Fermi-surface averages, Vph(q) => V p h, 
V c (q) V c J|, then there is only an s-wave solu- 
tion of Eq.([T]), A s , independent of k. The sound speed 
anisotropy actually changes the symmetry of the BCS 
state. While c/ is a constant in the isotropic medium, 
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FIG. 1: The Coulomb repulsion, /x m , as a function of the ratio 
of the electron wavelength to the screening length squared 
([3 = qg/2kp), and the electron- phonon coupling constant, 
A m , as a function of the ratio of the electron wavelength to 
the inter-plane distance squared, 7 = 7r 2 /2d 2 /c^(l + a) for 



a — 4 (inset) in s,p and d pairing channels. Here fi c 



it depends on the direction of q in any crystal. The 
anisotropy is particulary large in layered crystals like 
cuprate superconductors, where an elastic stiffness con- 
stant in the a — b plane is substantially greater than in 
the c direction (see [ilsl . [ItJ and references therein) . As 
an example, the measured velocity of longitudinal ultra- 
sonic waves along a — b plane, =4370 ms" 1 is almost 
twice larger than that along c axis, cj=2670 ins -1 in 
Bi 2 Sr 2 CaCu208+y It makes V p h(<l) anisotropic, 
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where a = (c 2 — Cj_)/cj_ is the anisotropy coefficient, 
which is about 2 in cuprates. Also the Coulomb repulsion 



is q dependent, V c (q) = Aire /Ve (q 2 



Here eg is the 



dielectric constant of the host lattice of the volume V , 
and q 2 = 8ire 2 N(0)/Veo is the inverse screening radius 
squared with the density of states (per spin), A^(0), at 
the Fermi surface. 

Now we can solve the master equation ([1} by ex- 
panding the order parameter and the potentials in sc- 
ries of the spherical harmonics. If the electron energy 
spectrum is quasi-two dimensional (2D) as in cuprates, 
£k ~ £ky , one can expand A(k) = ^2 m A m exp(im.</>) and 
Vph.ci^) = J2 m V phAl±i m ) exp[im(0 - </>')] in series of 
the eigenfunctions of the c-axis component of the orbital 
angular momentum, where <f> and </>' are the polar angles 
of k|| and kj. , respectively. 

The solution for the m-component of the order pa- 
rameter (to = 0, ±1,±2, ...) is found in the form A r 
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FIG. 2: The critical values of the electron-phonon coupling, A, 
for s-wave (gray color) and d-wave (a) and p-wave (b) Cooper 
pairs (bright colors), as functions of the parameters j3 and 
7. The Tolmachev-Morel-Anderson logarithm is set here at 
fi c \n(Lu p /Lur>) = 3, and a. = 3. 



ferent values of A„ and Am below and above the cut-off 
energy ll>d, respectively. Integrating in Eq.JT]) over £k', 
</>' and q± yields the following pair of equations, 
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Here A TO and fi m are the phonon-mediated attraction and 
the Coulomb pscudopotential in the TO-pairing channel, 
given respectively by 
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and 

Mm = Vj P dx[x + 13 + 1 - y/(x + I3){x + P + 2)} m 
Mc 2 J 

yjx(x + (3)(x + (3 + 2) 

(7) 

where A = N(0)C 2 /NMc 2 , 7 = n 2 /2d 2 k 2 F (l + a), 7 = 

7(1 + a), //c = 4e 2 d 2 iV(0)/7rl/eo, and /3 = <? 2 /2fc£, (note 
that A, He, and q s do not depend on the carrier density 
since iV(0) is roughly constant in the quasi-two dimen- 
sional Fermi gas). 

The effective attraction of two electrons in the Cooper 
pair with non-zero orbital momentum turns out finite at 
any finite anisotropy, a^O, but numerically smaller than 
in the s-channel, Fig.l (inset) , as is also seen from its 
analytical expressions for s-wave pairing, m — (A s ), p- 
wave pairing, m = 1 (Ap), and for d-wave pairing, m = 2 
(Ad), obtained by integrating in Eq.©. When the inter- 
plane distance is much larger than the wavelength of elec- 
trons, 7 <C 1, one obtains A s « A, A p Aa(7/2) 1/<2 /2, 
and Xd ~ A p . In the opposite limit, 7 » 1 one finds 
A s w A(l + a), A p « Aa(2/7) 1 / 2 /3, and X d w A p /5. 

The Coulomb repulsion turns out much smaller in the 
unconventional pairing states than in the conventional s- 
wave state, Fig.l, which is also seen from the analytical 
expression for /i m , Eq.([7]). If 7 < f3, the repulsion con- 
stant Mm drops as l/(3 m+1 in the m-channel at strong 
screening, when (3 I. It provides a wide region with 
unconventional pairs in the "7 - /3" parameter space, in 
spite of the lower values of their electron-phonon cou- 
pling constants, Fig.l. Indeed, the critical temperature 
for m-Cooper pairing is 

T cm = I.Ulud exp (- 1 ) , (8) 

V A m - Mm / 

where ii* n = At m /[l+/i m ln(wp/wi))], as found from Eq.Q 
and Eq.©. The m-pairing state with the lowest value 
of A//i c is the ground state of the BCS superconductor, 
where the critical ratio \ j \x c is determined by the condi- 
tion A m = Mm as the function of the parameters 7 and (3. 
The critical d-wave surface, \j \i c = <S(7, (3), defined by 
this condition, is found below s-wave and p-wave surfaces, 
if /5 > 1, so that the d-wave state is the ground state 
in this region, as seen from Fig. 2. Higher momentum 
states, m ^ 3, have even a smaller Coulomb repulsion at 
a large (3, Eq.([7]), so that they can be the ground state 
as well, if 7 is so small, that A m in Eq.([5]) is almost m- 
independent for m ^ 1. However, an in-plane anisotropy 
of the sound velocity, compatible with the symmetry of 
the perovskite lattice, could readily stabilise the d-wave 
state against higher-momentum states. Naturally, if the 
sound speed is enhanced along the diagonal directions of 
the in-plane primitive cell, the d-wave order parameter 
would be zero along diagonals of the 2D Brillouin zone, 
where it changes its sign. 

Using the simplest parabolic approximation for a 2D- 
electron energy spectrum we can draw some conclusions 
on the carrier-density evolution of the order-parameter 
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FIG. 3: Critical sound-speed anisotropy, q/(1 + q) 1//2 = 
(cy — c 2 L )/c^c±, as a function of doping, x, for A = 
(solid lines correspond to d and s states, and dashed line to 
p-state). With increasing carrier density there is a quantum 
phase transition at x — x c from a d-wave to an s-wave super- 
conductor, when a > a c , and two quantum phase transitions 
from d-wave to the normal state and from the normal state 
to the s-wave state when a < a c . 



symmetry. Within this approximation, k 2 F = 27rdn and 
AT(0) = mV/2i:dh 2 , where n = 2x/n is the carrier den- 
sity and x is the doping level as in La2- a; Sr 2 ;Cu04 with 
the unit cell volume f2. The ratio of the parameters 
(3 = me 2 n/2TTh 2 d 2 e x and 7 = irfl/8d 3 x « O.OU/x is 
independent of the carrier density, /3/j — 4me 2 d/ir 2 h 2 eo, 
which is approximately 5 for the values of m — 4m e 
and eq = 10. Fixing the value of the EPI constant 
at A — /i c /12 and taking /j, c \a(u)p/uiD) — 3, we draw 
the anisotropy- doping phase diagram, Fig. 3, with the 
critical lines for s, p and d order parameters, defined 
by A m = A*m- The state with the lowest magnitude of 
the anisotropy, a/(l + a) 1 / 2 , is the ground state of the 
BCS superconductor. At substantial doping the screen- 
ing length becomes larger than the typical wavelength of 
electrons, (3 — > 0, so that the s-wave state is the ground 
state at a large number of carriers per unit cell for any 
anisotropy. On the contrary, the Coulomb repulsion is 
reduced to the local interaction at a low doping, [3 — > 00, 
and d-wave Cooper pairs are the ground state even at 
very low value of the anisotropy, Fig. 3. Interestingly, 
s— and d-states turn out degenerate at some intermedi- 
ate value of doping, x = x c . Hence there is a quantum 
phase transition with increasing doping from d— to s- 
superconducting state, if a > a c , and from d— to the 
normal state and then to the s-wave superconductor, if 
a < a c , see Fig. 3. 

In the strong-coupling regime, A 1, the pairing is in- 
dividual [13|, in contrast with the collective Cooper pairs. 
While the Bose-condensate of individual bipolarons could 
break the symmetry on a discreet lattice [18|, [l9| , here I 
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propose a symmetry breaking mechanism, which works 
even in a continuum model, where the ground state, it 
would seem, be s-wave [Hj] to satisfy the theorem 

The unscreened Frohlich EPI with optical phonons in 
layered ionic lattices like cuprates significantly reduces 
the finite-range Coulomb repulsion, since eo > 1 , al- 
lowing the deformation potential, Eq.Q, to bind carri- 
ers into real-space bipolarons, if A ^ 1. The potential, 
V(r) — — J2q Vph(q) exp(iq • r) is non-local in real space, 



v(t) = -v ph n 



5M 
d 



4tt(1 



,)l/2 r 3 



(9) 



falling as 1/r 3 at the distance r ^> d between two car- 
riers in the plane, where V p h = C 2 /Mcf. While its lo- 
cal part (oc S(r) = S(x)S(yj) is negated by the strong 
on-site repulsion U, the second non-local part provides 
bound pairs of different symmetries with the binding en- 
ergies A s > A p > Ad > ... in agreement with the the- 
orem. However, there is the residual Coulomb repulsion 
between bipolarons, v c (R), significantly reduced by opti- 
cal phonons. If we approximate the bipolaron as a point 
charge 2e, then v c (R) w 4e 2 /eoi?. However since bipo- 
larons have a finite extension, £, there are corrections to 
the Coulomb law. The bipolaron has no dipole moment, 
hence the most important correction at large distances 
between two bipolarons, i? ^> £, comes from the charge- 
quadrupole interaction [l[, v c (R) — 4e 2 (l±r/£, 2 / R 2 ) / cqR, 
where r\ is a number of the order of 1, and plus/minus 
signs correspond to bipolarons in the same or different 
planes, respectively. The dielectric screening, eo is highly 
anisotropic in cuprates, where the in-plane dielectric con- 
stant, eo ii, is much larger then the out-of-plane one, eo± 
20]. Hence the inter-plane repulsion provides the major 



contribution to the condensation energy. Since £ 2 oc 1/A, 
the repulsion of unconventional bipolarons with smaller 
binding energies, Ad, A p < A s , is reduced compared with 
the repulsion of s-wave bipolarons. As a result, with in- 
creasing carrier density we anticipate a transition from 
BEC of s-wave bipolarons to BEC of more extended p— 
and d-wave real-space pairs in the strong-coupling limit. 

A great number of observations, in particular phase- 
sensitive experiments [2l|, point to the unconventional 
d-wave symmetry of cuprate su per conductors (for a re- 
view see [12] )■ Several authors [ill ] have remarked that 
superexchange, and not phonons is responsible for this 
symmetry breaking. Here I arrive at the opposite con- 
clusion. Indeed, superexchange interaction, J, is pro- 
portional to the electron hopping integral, t, squared di- 
vided by the on-site Coulomb repulsion (Hubbard U) , 
J = At 2 /U, estimated as J « 0.15 eV in cuprates [ill ]. 
This should be compared with the acoustic-phonon pair- 
ing interaction, V p h, which is roughly the Fermi energy, 
V p h « Ep « At in a metal [3], or the bandwidth squared 
divided by the ion — ion interaction energy of the order 
of the nearest-neighbour Coulomb repulsion, Mcf » V c 
in a doped insulator [151 ] . The small ratio of two interac- 
tions, J /Vph ~ t/U < 1, or J/Vph ~ V c /U < 1 and the 
giant sound-speed anisotropy [l|| [TtJ favor conventional 
EPI as the only origin of the unconventional pairing both 
in underdoped cuprates, where the pairing is individual 
[ill], and in overdoped samples apparently with Cooper 
pairs 0. 
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